clc;clear;clf;
p1=[];
s1=[];
for j=1:200
   for i=1:j
    p1(i)=(1/i^2)*sin(i)^1025;
   end
   sum(p1)
   s1(j)=sum(p1);
end
plot(s1)